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Abstract. Recent X-ray observations of microquasars and Seyfert galaxies reveal the broad emission lines in their 
spectra, which can arise in the innermost parts of accretion disks. Simulation indicates that at low inclination 
angle the line is registered by distant observer with characteristic two-peak profile. However, at high incl i nation 
angles ( > 85°) two additional peaks arise. This phenomenon was discovered by iMatt. Perola fc Stellal il993fl 
using the Schwarzschild black hole metric to analyze such effect. They assumed that the effect is applicable to 
a Kerr metric far beyond of a range of parameters that they exploited. We check and confirm their hypothesis 
about such structure of the spectral line shape for a Kerr metric case. We use no astrophysical assumptions about 
physical structure of the emission region except assumption that the region should be narrow enough. Positions 
and heights of these extra peaks drastically depends on both the radial coordinate of the emitting region (circular 
hot spot) and the inclination angle. We find that best conditions to observe this effect are realized at ^ > 85° 
and r > 5rg and may exist in microquasars or low-mass black holes in X-ray binary systems, because there is 
some precession (and nutation of accretion disks) with not very long time periods (see, for example, SS433 binary 
system). The line profiles for different inclination angles and radial coordinates are presented. To analyze an 
influence of disk models on the spectral line shapes we simulate the line profiles for Shakura - Sunyaev disk model 
for accretion disks with the high inclination. 
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1. Introduction 

More than ten years ago it was predicted that pro- 
files of lines emitted by AGNs and X-ray binary 
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erences therein). Recent X-ray 
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(11993D). Generation of the broad fluorescence 
lines as a result of irradiation of cold accretion 
disk was discussed by many a ut hors (see, fo r exam - 
ple. iMatt. Perola fc Pirol lll99ll): iMatt et aJ lll992la[) : 
iMatt fc Perolal l)l992() : iMatt. Fabian fc Ros^ l|l993|) : 



1120021) and references therein) confirm these consider- 
ations in general and reveal broad emission lines in 
their spectra with cha racteristic two-peak profiles. A 
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^ So me of th e m are microquasars ( f or details see, for ex- 
ample , iCreined |l99i); iMirabell J2000l) : iMirabel fc Rodrigued 
J2002D '). 



comprehensive review bv lFabian et al.l 1 20001) summarizes 
the detailed discussion of theoretical aspects of possible 
scenarios for generation of broad iron lines in AGNs. 
These lines are assumed to arise in the innermost parts of 
the accretion disk, where the effects of General Relativity 
(GR) must be taken into account, otherwise it appears 
very difficult, if any, to find the natural explanation of 
the observed line profile. 
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Numerical simulations of the line structur e coul d 
be f ou nd in a number of p ap ers: llKqiimal 1 1991 k 



LaoJ lll99lh: iBao & Stuchlik I lll992l): iBad Jl99i: 



Bao. Hadrava fc OstgaardI ll994l): iRauch fc BlandfordI 



_,tg^ 

19931 11994): iBromlev. Chen fc Milled 
^antonetalJ il9 9/l: 'Paricv fc BromlevI lll£ 




Pariev. Bromlev fc Miller. (,2001.) : iRuszkowski MJ pOOQl : 
mJ |2002|) '). Thev indicate that the accretion disks in 
Seyfert galaxies are usually observed at the inclination 
angle 9 close to 30° or less. It occurs because according 
to the Seyfert galaxy models, the opaque dusty torque, 
surrounding the accretion disk, so, such structure does 
not allow us to observe the disk at larger inclination 
angles. 

However, at large inclination angles 9 > 80° 
the new observational manifestations of GR could 
arise, (^att. Pcrola fc Stella (1993) discovered such phe- 
nomenon for a Schwarzschild black hole, moreover the au- 
thors predicted that their results could be applicable for a 
Kerr black hole over the range of parameters which were 
exploited). The authors mentioned that this problem was 
not analyzed in details for a Kerr metric case and it would 
be necessary to investigate this case. Below we do not use 
some specific model on surface emissivity of accretion (we 
only assume that the emitting region is narrow enough). 
But general statements (which will be described below) 
could be generalized on a wide disk case without any prob- 
lem. Therefore, in this paper we check and confirm their 
hypothesis for Kerr metric case and for a Schwarzschild 
black hole using another assumptions about surface emis- 
sivity of accretion disk. In principle, such phenomenon 
could be observed in microquasars and X-ray binary sys- 
tems where there are neutron stars and black holes with 
stellar masses. 

In Subsection 3° we discuss disk models used for simu- 
lations. In Subsection 4° we present the results of simula- 
tions with two extra peaks in the line profile. These extra 
peaks exist because of the gravitational lens effect in the 
strong field approximation of GR. In Subsection 5° we dis- 
cuss results of calculations and present some conclusions. 



2. Numerical methods 

We used a n approach which was discussed in details in 
papers hv I Zakharovl ^19911 11991: Fz'akhaTOv fc T^eninl 
lll999[ '200023^: 'Zakh a,rov fc Remn I ll2nn2bl ^QQ^; 
IZakharov et al, (,,2002') . The approach was used in par- 
ticular to simu l ate s pectral line shapes. For example, 
IZakharov et al.l ((20021 ) used this approach to simulate an 
influence of magnetic field on spectral line profiles. This 
approach is b ased on results of qualita tive analysis (which 
was done by I Zakharovl (|l98l ll98i) for different types 
of geodesies near a Kerr black hole). The equations of 
photon motion in Kerr metric are reduced to the following 
system of ordinary differential equations in dimcnsionlcss 
Boyer - Lindquist coordinates (|i Zakharov. (,199L .1994) : 
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where A = - 2r a^; = Q/M'^E'^ and £ = L JME 
are the Chandrasekhar constants (|cLandrasekhi3 dl 9831) 1 
which are derived from the initial conditions of the emitted 
quantum in the disk plane. The system (E3)-(ini) has two 
first integrals 
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which can be used for the accuracy control of computation. 

Solving Eqs. iQl"® for monochromatic quanta emit- 
ted by a ring we can calculate a specral line shape I^, (r, 9) 
which is registered by a distant observer at inclination 
angle 9. 

The numerical integration has been performed us- 
ing_the combination of the Gear and Adams methods 
dGead (Il97ll)') and reahzed as the standard package by 
iHin drnarshl (Il983l) : IPetzoldl lll983l) : iHiebert fc Shamping 
(198(3 ). We obtain the entire disk spectrum by summation 
of sharp ring spectra. 

3. Disk model 

To simulate the structure of the emitted line it is nec- 
essary first to choose a model for an emissivity of an 
accretion disk. We exploit two different models, namely 
we consider narrow and thin disk moving in the equa- 
torial plane near a Kerr black hole as the first model 
and as we analyze the inner wide part of accretion disk 
with a tem perature distributioii which is chosen accord- 
ing to the I Shakura fc SunvaevI (figySl with fixed inner 
and outer radii and Tq as the second model. Usually 
a power law is us e d for a wide disk emissivi t y (see , 
for example. iLaoil (Il99lh- iMatt. Perola fc Pirol 'lll99lh: 
IMartoccchia fc MattI lll99f)l): IMartoccchia. Karas fc Mat!il 



(|200(J) : IM artoccchia. Matt fc Karas I (|2002l) ) . However, 
another models for emissivity could not be forbidden for 
so wide class of accreting black holes, therefore, just to 
demonstrate how another emissivity law could change line 
profiles we exploit such emissivity law. 
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So, at the first, we assume that the source of the emit- 
ting quanta is a narrow and thin disk rotating in the equa- 
torial plane of a Kerr black hole. We also assume that the 
disk is opaque for radiation, so that a distant observer sit- 
uated on one disk side cannot register the quanta emitted 
from its other side. 

For the sake of computational simplicity we suggest 
that the spectral line is monochromatic in its co-moving 
frame. It can really be adopted because even at T = 10® K 
the thermal line width 




appears to be much less than the Doppler line width as- 
sociated with the Kepler velocity of the disk rotation. 

For second case of the Shakura - Sunyaev disk model, 
we assume that a local emissivity is proportional to surface 
element and and T^, where T is a local temperature. The 
emission intensity of the ring is proportional to its area. 
The area of emitting ring (width dr) differs in the Kerr 
metric from its classical expression dS = 2i:rdr and should 
be replaced with 



27r (r^ +a'^) 
dS = ^ -!— dr. 



(9) 



Thus, the total flux density emitted by the disk and 
registered by distant observer is proportional to the inte- 
gral 



MO)^ / lAr,0)T\r)dS, 



(10) 



where (r, 9) is obtained from the solution of equations 
(C3)-lini), T{r) - from the appropriate dependence for a hot 
disk and dS - from Eq. ®. 

For simulation we assume that the emitting region lies 
entirely in the innermost region of a-disk (zone a) from 
fout — 10 to Tin — Svg and the emission is monochro- 
matic in the co-moving frame. ^ The frequency of this emis- 
sion set as a unity by convention. 

4. Simulation results 

Spectral line profiles of a narrow ring observed at large in- 
clination angles 9 > 85° and different radii r are shown in 
Fig.n(one could say that is a Kerr metric generalization of 
Fig. 1 from lMatt. Perola fc Stellal l)l993|) which was drawn 
using calculations for the Schwarzschild case) . The ring is 
assumed to move in the equatorial plane of a Kerr black 
hole with almost extreme rotation parameter a = 0.9981. 
The inclination angle increases there from left to right and 
radial coordinate - from bottom to top. The figure indi- 
cates that practically there are no new specific features of 
profiles, thus, the line profile remains one-peaked with a 
maximum close to 1.6 Eiat and very long red wing without 

^ We use as usual the notation Vg = 2GM/c^. 



any significant details. For lowest radii there are no sig- 
natures of multiple peaks of spectral line shapes even for 
high inclination angles (the bottom raw in Fig. 1 which 
corresponds to r — 0.8 r^). 

With increasing the radius to r = 1.2rg the additional 
blue peak arises in the vicinity of the blue maximum at the 
highest inclination angle 9 = 89° . The red maximum is so 
small that no details can be distinguished in its structure. 
At lower inclination angles 9 < 88° the blue maximum also 
has no details and the entire line profile remains essentially 
one-peaked. 

For r = 3 the additional detail in the blue peak 
appears for 9 > 85°. Thus, for 9 = 85° we have a fairly 
distinguished bump, for 9 = 88° it changes (transforms) 
into a small complementary maximum and for 9 = 89° 
this maximum becomes well-distinguished. Its position in 
the last case differs significantly from the main maximum: 

E3 = 1.12 Eiab, E4 = 1.34: Eiab- 

With further increasing the radius the red maximum 
bifurcates also. This effect becomes apparently visible for 
r = 5rg and 61 > 88°. Thus, for 9 = 85° we have only a 
faintly discernible (feebly marked) bump, but for 9 = 88° 
both complementary maxima (red and blue) arise. For 9 = 
89° we have already four maxima in the line profile: Ei — 

0.63 Eiab, E2 = 0.66 Eiab, E3 = 1.07 Eiab, Ei = 1.28 ^,afc. 

Note that the splitting for the blue and red maxima is not 
equal, moreover, E4 — E3 ^ 7 {E2 — Ei). 

For r = 10 Tg the profile becomes more narrow, but the 
complementary peaks appear very distinctive. We have 
four-peak structure for 9 > 88°. It is interesting to note 
that for 9 — 89° the energy of the blue complementary 
peak is close to its laboratory value. 

Note that the effect almost disappears when the radial 
coordinate becomes less than Vg, i.e. for the orbits which 
could exist only near a Kerr black hole. 

Fig. 121 demonstrates the details of the spectrum pre- 
sented in the top row in Fig. Thus, the left panel in- 
cludes all the quanta emitted by a hot spot at r = 10 rg 
with 9 > 89° at infinity (the mean value could be counted 
as 85? 5 if the quanta distribution would be uniform there) 
but with much higher resolution than in Fig. ^ The 
spectrum has four narrow distinctive maxima separated 
by lower emission intervals. In the right panel which in- 
cludes all the quanta with 9 > 85? 5 the right complemen- 
tary maxima is even higher than the main one. The blue 
complementary maxima remains still lower than its main 
counterpart, but it increases rapidly its intensity with in- 
creasing the inclination angle. It follows immediately from 
the comparison of left and right panels. Some "oscillation 
behavior" of the line profile between the maxima has a 
pure statistical origin and does not refer to physics. 

As an illustration a spectrum of an entire accre- 
tion disk at high inclination angles in the Schwarzschild 
metric is shown in Fig. |31 (In reality we've calculated 
geodesies for the quanta trajectories in the Schwarzschild 
metric using the same Eqs. as for a Kerr met- 

ric, but assuming there a = 0.01.) The emitting re- 
gion (from 3 to 10 rg) lies as a whole in the in- 
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Fig. 1. Line profiles with specific details, which arise for high inclination angles 6 > 85° due gravitational lens effect in 
the strong gravitational field approach. The Kerr (rotation) parameter was chosen like a = 0.9981. The radii decrease 
from top to bottom, the inclination angles increase from left to right; their values are shown in each panel along with 
the number of quanta included in spectrum. 




Fig. 2. Details of a hot spot line profile for the most distinctive case with r = lOvg and a — 0.9981 (see the top row 
of Fig.^l. The images of all orders are counted. Left panel includes all the quanta, registered at infinity with 9 > 89°. 
Right panel includes the quanta with 9 > 89? 5. 




Energy Energy 



Fig. 3. Details of a line structure for a-disk in the Schwarzschild metric with outer and inner edges of emitting region 
equal to rout = lO'^g and ri„ — Srg, respectively. Left panel includes all the quanta, registered at infinity with 9 > 89°. 
Right panel includes the quanta with 9 > 89? 5. 



descrip t ion o f 

this model was given by I Shakura fc SunvaevI l|l973() : 



nermost region of a-disk (the detailed 
y 

iLipunova fc Shakural l)2002|] 'l. 



As it follows from the Fig. the blue peak may con- 
sists of the two components, whereas the red one remains 
unresolved. 



5. Discussion and conclusions 



The complicated structure of the line profile at large incli- 
nation angles is explained by the multiple images of some 
pieces of the hot ring. One could point out that the result 
was obtain in the framework of GR without any extra 
physical and astrophysical assumptions about character 
of radiation etc. For a Kerr black hole we assume only 
that radiating ring is circular and narrow. 
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The problem of multiple images in the accretion 
disks and extra peak s has first beer i considered by 
I Matt. Perola fc Stellal lll993tl see also iBao fc Stuchlik 
I 1992|) : iBad ijlQQ^ iBao. Hadrava fc Ostgaarc 
( 1994i 1 199611). Us ing numerical simulations 
iMatt. Perola fc Stellal l)l993|) proved the statement 
for the Schwarzschild metric and suggested that the 
phenomenon is applicable to a Kerr metric over the 
range of parameters that the authors have analyzed. 
iMatt. Perola fc Stellal l|l993(l noted also that it is nec- 
essary to realize detailed calculations to confirm their 
hypothesis. We verified and confirmed their conjecture 
without any assumptions about a specific distribution of 
surface emissivity or accretion disk model (see, Fig. lll2|l . 

We confirmed also a conclusion of 
iMatt. Perola fc Stellal l(l99.'^ that extra peaks are 
generated by photons which are emitted by the far 
side of the disk, therefore we have a manifestation of 
gravitational lensing in the strong gravitational field 
approach for GR. 

Some possibilities to observe considered fea- 

tures of spectral U n e pro files wa s con sidered by 

iMatt. Perola fc Stellal lll99.1) : iBaol lll99.ll) . The au- 
thors argued that there are non-negligible chances to 
observe such phenomenon i n some AG Ns and X-ray 
binary systems. For example, iBaol l)l993(l suggested that 
NGC 6814 could b e a candidate t o deni onstrate such 
phenomenon (but iMadeiskT et all l)l993|) found that 
the peculiar properties of NGC6814 are caused by a 
cataclysmic variable like an AM Herculis System). 

However, it is clear that in general a probability to ob- 
serve objects where one could find such features of spectral 
line is small. Moreover, even if the inclination angle is very 

o 

close to 90 , a thickness of disk ( shield or a torus around 
accretion disk) may not give us a possibility to look at the 
inner part of accretion disk. Let us discuss astrophysical 
situations when the the inclination angle of accretion disk 
is high enough. 

About 1% of all AGN or microquasar systems 
could have an inclination angle of accretion disk > 
89 . For example, fKormcn dv et alJ l)l996ll found that 

o 

NGC3115 has very high inc lination angle about 81 
i Kormendv fc Richston'el lll992^ discovered a massive dark 
object M ~ 10^ M0 (probably, massive black hole) in 
NGC3115). Perhaps, we have much higher probability to 
observe such phenomenon in X-ray binary systems where 
black holes with stellar masses could be. Taking into ac- 
count an precession which actually observed for some X- 
ray binary systems (for example, there is a significant 
precession of the accr etion disk for SS433 binary system 
IChereoashchukl l)2002l) .'^ moreover since the inclination of 

o 

the orbital plane is high (i 79 for this object), so, some- 
times we could have a possibility to observe almost edge- 
on accretion disks for such objects. Observations indicated 

^ f Shakura ( 197i) predicted that if the plane of an accretion 
disk is tilted relatively to the orbital plane of a binary system 
the disk can precess. 



that there is a strong evidence that optically bright accre- 
tion disk in SS433 is in supercritical regime of accretion. 
Fir st description for supe r critica l accretion disk was given 
by I Shakura fc SunvaevI l)l973(l . Even now such model 
is discussed in to exp lain observational data for SS433 
IChereoashchukl l)2002j) . We used a t emperature distribu- 
tion fr o m Shakura - Sunyaev model \ Shakura fc SunvaevI 
(197^; iLipunova fc Shakural l|2002l) ) for inner part of ac- 
cretion disk to simulate shapes of lines which could be 
emitted from this region (Fig. OJ. Therefore, we should 
conclud e that the properties of s pectra l line shapes discov- 
ered bv iMatt. Perola fc Stellal 11993) are confirmed also 
for such emissivity (temperature) distributions which cor- 
respond to Shakura - Sunyaev model. 

Thus, such properties of spectral line shapes are robust 
enough in respect to wide variations of rotational param- 
eters of black holes and a surface emissivity of accretio n 
disks as it was predicted bv lMatt. Perola fc Stellal 119931) . 
In this work their conjecture it was verified and confirmed 
not only for a Kerr black hole case but also for another 
dependence of surface emissivity of accretion disk. 
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